Spin ordering in magnetic quantum dots: Prom core-halo to Wigner molecules 



R. Oszwaldowski, 1 P. Stano, 2 A. G. Petukhov, 3 and Igor Zutic 1 

1 Department of Physics, University at Buffalo, Buffalo, NY 14-260-1500 
2 Institute of Physics, Slovak Academy of Sciences, 845 11 Bratislava, Slovakia and 
Department of Physics, University of Basel, Klingelberstrasse 82, 4056 Basel, Switzerland 



o 

(N 
O 

Q 



c3 



-a 
c 
o 

(N 
> 

(N 



(N 



X 



Department of Physics, South Dakota School of Mines and Technology, 



City, SD 57701 



The interplay of confinement and Coulomb interactions in quantum dots can lead to strongly 
correlated phases differing qualitatively from the Fermi liquid behavior. We explore how the presence 
of magnetic impurities in quantum dots can provide additional opportunities to study correlation 
effects and the resulting ordering in carrier and impurity spin. By employing exact diagonalization 
we reveal that seemingly simple two-carrier quantum dots lead to a rich phase diagram. We propose 
experiments to verify our predictions, in particular we discuss interband optical transitions as a 
function of temperature and magnetic field. 
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With high tunability of their parameters, quantum 
dots (QDs) are ideal systems for exploring correlation 
effects^— While in 3 dimensions the correlation-induced 
Wigner crystal^ is elusive and only expected in the limit 
of an extremely low carrier-density its nanoscale ana- 
log, the Wigner molecule (WM),-^ was observed in QDs 
at much higher densities*^ An increase in the relative 
strength of Coulomb interactions qualitatively changes 
the liquid-like independent-particle picture to that of 
WM characterized by electron localization and strong an- 
gular order ,2iL& 

We expect that doping QDs with magnetic impuri- 
ties^— (typically Mn) opens unexplored opportunities 
to study the nanoscale correlations. Through Mn-carrier 
exchange interaction, the correlations can be enhanced, 
imprinted on Mn-spins, and thus observed. Several key 
advances in elucidating correlations in non-magnetic QDs 
were accomplished in two-electron systems i 2 -^^ How- 
ever, even in simple circular QDs, identifying WMs is 
complicated by insufficient accuracy in the treatment 
of corrclation o 2 ' 14 and artifacts of the mean-field and 
Hartree-Fock approaches^ To understand these sys- 
tems, the exact diagonalizatio n 8 ! 9 ' 14 is particularly suit- 
able, corroborating analytical findings for two-electron 
correlations in QDsj 16 ' 17 

Here, we generalize the exact approac h 18 ' 19 to probe 
the charge and spin densities of carriers as well as 
Mn-spin ordering in magnetic QDs with two carriers 
(holes) j2£ To elucidate the stability of the magnetic or- 
dering, we consider different Mn-doping configurations, 
deformation of circular QD confinement, and examine 
the effects of temperature T and magnetic field B. 

The phase diagram in Fig. [TJ for Mn-doped circu- 
lar QDs, shows three magnetic groundstates, to be con- 
trasted with the spin-singlet groundstatc of non-magnetic 
QDsi£ The groundstate phase changes with the doping 
radius i? e ff and the fraction, iMn, of cations replaced 
by Mn-atoms in the QD. The carrier spin density is im- 
printed on Mn-spins forming three patterns correspond- 
ing to: two pseudo-singlets (PS) characterized by total 
hole-spin zero, but non-zero hole-spin density^ and a 



spin triplet (T), see Fig. HJa-c) 
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FIG. 1. (color online) Groundstate phase diagram (PS ver- 
tical and skew, T horizontal hatching) as a function of Mn- 
content xmyi and doping radius R e g, for a double-occupied 
circular dot, Huj m , y — 25 meV, at zero temperature. Insets 
(a-c): QD-top view (in-plane coordinates in nm) of Mn-spin 
patterns; spins up: light, down: dark, (a) Core-halo (CH): 
RsB = 2 nm, im u = 2%, (b) triplet (T): R c g = 3.5 nm, 
iMn = 3%, (c) spin- Wigner Molecule (sWM): R c g = 3.5 nm, 
iMn = 1%, (d) Hole-spin density for the pattern in (a).— 



Before we provide a detailed analysis, it is instructive 
to view the PS state in Fig.QJc) as the spin-WM. In non- 
magnetic circular QDs, their WM "dimer-like" phased 
can be fully revealed only in the pair-correlation func- 
tion*^ However, a similar phase can be directly detected 
in magnetic QDs: the Mn-pattcrn of spin-WM reflects a 
double-peaked hole-spin density. The separation between 
carriers, characteristic for WMs^ is enhanced with Mn, 
which provides a spin structure. 

We use the total QD Hamiltonian, H = Hj + H cx , with 
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typical 2D non-magnetic and exchange parts^i where 
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(i) 



the holes are labeled by i, m* is the effective mass, e the 
electron charge, and e the dielectric constant, momentum 
7Tj includes the vector potential of field The p — 

d exchange interaction between Mn-spins and confined 
holes has the Ising form 24 
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because of the strong z-axis anisotropy, arising from spin- 
orbit interaction in the 2D QDs. Here, f3 is the exchange 
constant, s z is the heavy-hole pseudospin with projec- 
tions s z = ±3/2, while S z are operators of z-projection 
of Mn-spins positioned at Rj, and N is the number of 
Mn-spins in the dot. 

Since H cx does not contain spin-flip processes, the 
total wavefunction is a product of the hole and Mn- 
spin parts. The partition function of the system can 
be calculated using Gibbs canonical distribution, Z = 
Tt s jz E„ cx P [~ E n ({Sjz}) /k B T] , where k B is the Boltz- 
mann constant, and the hole eigenvalues, E n , depend on 
N numbers Sj Z , with each Sj Z = — S, S = 5/2 (index n 
runs over hole states for fixed {Sj Z }). To calculate Z, one 
would need to solve Q N replicas of the hole Schrodinger 
equation, with N ~ f 2 — fO 3 . 

We can overcome this obstacle of computational com- 
plexity by partitioning the dot into N c square cells, each 
containing few Mn-spins, and neglecting spatial variation 
of the two-hole wavefunction <£>„ within each cell through 
the use of the average hole-spin density 



1 



h z N k 



jeN k 



^2s iz 6(r l - Rj) 



(3) 



where h z is the QD height. For a given cell with Nk 
spins creating a magnetic moment M k , a distribu- 
tion function of the average dimcnsionlcss magnetization, 
m k = — A/fc/gMnMs-^fc, can be expressed as Y{mk) oc 
cxp [—Gk{rrik/S)/k B T]. Here jMn = 2 is the Mn g-factor. 
The Gibbs free energy of the Nk non-interacting spins, 
Gk{mkl 'S), is obtained by Legcndre transformation^ 



G k (x) 
N k k B T 



sinh [(l + l/2S)Bg\x)] 



sinh [(l/2S)Bg\x)] 
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where Bg is the inverse of the Brillouin function B$. 
Using Y(rrik), we transform Z with exponential accuracy 
as 
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cxp[-Gl l ot /k B T]l[dmk, (5) 



fe=l 



where G^ot = E fe G k (m k /S) + E n ({m k }). For 
any n, each integral in Eq. ([5]) can be evaluated 
using the steepest descent method. Equation for 
the saddle point, combined with the Hellmann- 
Feynman theorem, -/3N k (s k )/3 + SMn-^fc/XB-B = 
dEn(mi,...,m k ,...,m,N c )/dmk, leads to the self- 
consistency condition 

m fc = SB S [S (p (s k ) n /3 - g M nLi B B) /k B T] . (6) 

Our analysis shows that the self-consistency condition 
Eq. ^ depends on the quantum-mechanical average, 
(sk) n , relevant for small systems such as QDs, rather 
than on the thermal average,— thereby avoiding arti- 
facts arising from imposing the thermodynamic limit on 
a nanoscale system. We use Eq. (O to find a global min- 
imum of G^ ot ({m k }) corresponding to the groundstate, 
where stands for PS or T in different regions of the 
phase diagram. 

We find the eigenstates of H self-consistently: For fixed 
values of m k (randomly initialized), we obtain interme- 
diate two-hole states. Since [H, E 2 ]=0, where S z is the 
total hole-spin z-projection, the states are either in PS 
(£ z = 0) or T (£ z = ±3) orthogonal subspaces^i We 
choose the lowest state in each subspace, use Eq. (|6]) to 
obtain new m^ , and restart exact diagonalization, iterat- 
ing until convergence. The groundstate nature, either T 
or PS, depends on the QD parameters (for ime = = B 
the groundstate is a singlet). 

We use ZnTe parameters: hole mass m* = 0.2 electron 
mass^ e = 9.4eor^ eo the vacuum permittivity, Nq(3 = 
— 1.05 eVr^ and N = A/a 3 the cation density with the 
lattice constant a ~ 6.1 A. For typical self-assembled 
QDs: hiij x ^y = 10 — 30 meV, im e < 5%, our standard 
values are iMn = 1%, h z = 1.8 nm, Huiq = hoj x , y — 
25 meV, corresponding to a characteristic length 2 ^ Iq = 
(h/m*uJo) 1/2 = 3.9 nm. 

We now revisit Mn-patterns in Fig. [TJ calcu- 
lated at T = 0, with a doping profile iMn{l + 
exp [(r — -R ff) /?]} _1 of width £ = 0.25 nm, and radius 
i? e fi , here r is the distance from the center. In addition 
to the triplet groundstate (all Mn-spins parallel) at rela- 
tively large R a f{ and xme, and the spin-WM groundstate 
at smaller iMn, another PS forms, Fig. [IJa), we term it 
"core-halo" i2S Inset (d) shows the resulting average spin- 
density (sfc). In contrast to spin-WM, the Mn-pattern 
and (sfc) for core- halo preserve the circular symmetry of 
the QD. For £ — > (no Mn for r > R e g) and small enough 
R e s (e.g. -R e ff = 2 nm, XMn = 1%) the core-halo patterns 
become purely "ferromagnetic." We find that such uni- 
directional patterns form also for other inhomogeneous 
^Mn-profiles, confirming predictions from Rcf.i^i 

In the i? e ff — > oo regime, we study stability of T and 
PS, termed magnetic bipolarons* 2 ^ Figure [D(a) shows 
that, already at modest iMn, the Mn-induced energy gain 
of PS exceeds that of a magnetic polaron forming for a 
single carrier^ This suggests that the robustness of mag- 
netic bipolarons is similar to that of the well-established 
single magnetic polarons. 
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FIG. 2. Temperature evolution of pseudosinglet Eps (green) 
and triplet Et (red) energies for a circular QD, hu) XtV = 25 
meV, and iMn = 1%. Bold tics show the high-T (zero p — d 
exchange) limit, (a): Mn-induced energy gains at T = 
K. Dots/solid line: numerical/variational results^ for PS, 
A_Bps = Eps (imb = 0) — Eps ■ Dash-dotted line: single mag- 
netic polaron A£mp = iMn^o \/3\ 5/2. A£W < AEps for 
iMn > 0.7%. Insets (b-d): rrik in units of 1/4 as a function of 
position (nm); (b) triplet at 1 K, (c,d) PS at 1 K, 3 K. 



Figure [5] shows the temperature dependence of bipo- 
laron energies. Owing to the stronger exchange interac- 
tion, the triplet approaches the high-T limit at a higher 
T than PS. The asymptotic trend — 1/T in E^r, typical 
for paramagnets, is expected since the effective exchange 
field of triplet holes, acting on Mn, is nearly independent 
of the Mn-spin alignment. In contrast, PS has a second- 
order transition to singlet^ Insets (b-d) show the mag- 
netization [Eq. (j5])] of PS and T . Finite temperature has 
a different effect on the two states. For PS, nik retains its 
overall K shape, and decreases uniformly. In contrast, 
the saturated Mn-pattern of T, Fig. QJb), undergoes a 
transition to a state with spontaneously broken circular 
symmetry: two symmetrical and unidirectional peaks ap- 
pear, Fig. &b). The peaks reflect the hole-spin density, 
which maximizes the Mn-induced energy gain through a 
linear combination of the two triplets with opposite an- 
gular momenta^ 

Having established the presence of magnetic bipolarons 
and Mn-patterns in circular QDs, it is crucial to ex- 
amine the more realistic case of asymmetric confine- 
ment. We introduce in-plane asymmetry through w x ,y = 
H/(m*l% y ), where l x>y = lo (1 + d) 1 M- The asymmetry 
stabilizes the spin-WM along the weaker-potential axis, 
Fig. [3ja), since separating the holes in the direction of 
the "softer" potential costs less potential energy. As in 
circular QDs [Fig. [TJ large R e s), the triplet becomes the 
groundstate with increasing iMn- For increasing asym- 
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FIG. 3. (color online) Groundstate phase diagram (PS verti- 
cal, T horizontal hatching) as a function of confinement asym- 
metry d and iMn, at T = 0. Insets: hole spin density of PS 
(a) and T (b), huo = 25 meV, d = 1, and XMn — 1%. 



metry d, the non-magnetic singlet-triplet gap decreases, 
so that T becomes the groundstate at lower iMn^ 

The results so far can be understood introducing an 
effective, spin-dependent interaction between the two 
holes. Each hole tends to polarize Mn-spins within an 
area defined by temperature and the confinement. Con- 
sequences of the overlap of the areas depend on relative 
hole spins. If they are antiparallel (PS), mutual repul- 
sion arises^ which would lead to complete separation 
of holes in absence of the confinement (magnetic bipo- 
laron would not form). In contrast, holes with parallel 
spins (T) effectively attract each other,— as they bene- 
fit from sharing the cloud of polarized Mn-spins each of 
them carries^ 

Opposite limits of the doping radius i? e ff offer a simple 
insight as to which of the two PS patterns is the ground- 
state. For R c h — > 0, the p — d exchange energy gain for 
spin-WM is zero, since this state produces (su) = at 
the QD center. On the other hand, for R c h — > oo and 
for weak confinement, the proper limit is that of two sep- 
arate, localized magnetic polarons, a scenario consistent 
with spin-WM rather than core-halo^ 

To corroborate our predictions, as well as other 
works considering magnetic interactions in closed-shell 
QDs ) 13 ' 21 i 36 ' 37 we propose experiments that test the ex- 
istence of magnetic bipolarons, and discriminate the dif- 
ferent Mn-patterns. One such probe is interband pho- 
toluminescencei 10 ' 11 ^ 8 With a sufficiently intense excita- 
tion, two kinds of emission lines appear, corresponding 
to 2 — >• 1 and 1 — >• QD-occupancy transitions. We 
calculate the photon energies, i? p h, for the standard pa- 
rameters, assuming: (i) type-II band alignment ) 28 i 32 (ii) 
Mn-spin pattern does not change during a recombina- 
tion event, (iii) the system recombines from its (two- or 
single-hole) groundstate. We show that E p h dependence 
on T and B allows to identify the bipolarons. 

We first consider varying T, Fig. Single polaron, 
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characterized by a unidirectional Mn-pattcrn, shows a 
1/T rcdshift. In contrast, thermal disruption of PS oc- 
curs with an abrupt change of the slope of Ep^ at few 
kelvin. 
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FIG. 4. (color online) Temperature dependence of the PS — > 1 
hole (solid), and 1 — > (dash-dotted) transition energies -E p h, 
for B = 0. To better compare the dependencies, each line is 
shifted: AE ph = E ph (T) - E ph (T -> oo). 

We next consider B \\ z (Faraday configuratio n 28 ' 39 ) . 
For a triplet, rrik \\ —B everywhere. The small change of 
A-Eph, (Fig. [SJ dashed) is mainly due to orbital effects, 
since mu , highly saturated at the low-T, is not very sen- 
sitive to B. For PS, a part of Mn-spins is aligned unfa- 
vorably, i.e., parallel to B. Increasing B changes their 
projection [Fig. (5][a,b)]. The accompanying change of 
Eps (not shown) becomes abrupt at a threshold B (the 
Mn-pattern becomes unidirectional at Bq), and then flat- 
tens out close to its asymptotic value: the non-magnetic 
singlet energy 42, For B > B Q (B ~ 1.5 T in Fig. [5J, 
the independence of PS — > 1 transition energies is dom- 
inated by changes of final-state energies, which react to 
increasing saturation of rrik- The line is split by the field, 
as the remaining hole can end up in two opposite spin 
states. No such splitting occurs for the T — > 1 and 1 — > 
transitions. Thus, the splitting, and the abrupt 'melting' 
at Bq, signal the PS groundstatc. 

Finally, the two possible PS states can be resolved us- 
ing selection rules for photoluminescence4i The circu- 
larly symmetric core-halo pattern forbids recombination 
with p-like excited electron states. However, such transi- 
tions are (weakly) allowed for spin-WM^ and would ap- 
pear at a strong optical pumping J 2 - As an alternative to 
photoluminescence, scanning tips with NV-centers could 
offer sufficient spatial sensitivity to probe the Mn-spin 
patterns 4 s - 



We expect that our findings will motivate future efforts 
to probe magnetization patterns and correlation effects 
at the nanoscale. Previously unexplored regimes are af- 
forded with Mn-doping. Effective internal magnetic fields 
in colloidal QDs can reach ~ 100 T4i beyond what is 
feasible with applied static magnetic fields. Consistent 
with recent advances in the field of nanomagnctism4^ an 
increasing number of experimental probes are likely to 
meet the challenge of detecting the predicted spin order- 
ing even in single QDs. 

Nuclear spins in few-electron III-V QDs could pro- 
vide a magnetically-active feedback similar to that stud- 
6p i ■ ■ -3-T 
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FIG. 5. (color online) B\\z dependence of the PS — > 1 (solid 
lines) and T ^ 1 (dashed), for T = 1 K. To avoid PS-T 
crossing, we set xyin = 0.5%. Each line is shifted: A_E p h = 
E ph (B) - E ph (B = 0). Insets: PS Mn-patterns for 0.4 (a), 
and 1.4 tesla (b); light (dark) for rrik > (rrik < 0). 



ied here. While the electron-nuclear spin interaction is 
weak, leading to a much smaller temperature scale for 
analogous polaron objects, this scale is known to be 
very strongly enhanced by electron-electron interactions 
in low-dimensional systems.— 

This work was supported by DOE-BES, ONR, meta- 
QUTE ITMS NFP 26240120022, CE SAS QUTE, EU 
Project Q-essence, APVV-0646-10 and SCIEX. 
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Spin Ordering in Magnetic Quantum Dots: From Core-halo to Wigner Molecules 

(Auxiliary Material) 



This material provides details helping to explain some statements in the main text. The 1st occurrence of the 
EPAPS reference in the main text points to Sec. IA11 the 2nd and 4th to Sec. Al.l, the 3rd to Sec. A1.2, the 5th and 
6th to Sec. H2J 



Al. VARIATIONAL CALCULATION 

We build the pseudo-singlet (PS) for the spin-WM 
ground state as 

$ps = ^ [fu{ri) tpd(r 2 ) xt (1) (8 Xi (2) 
-Puto) <Pd(ri)xt ( 2 ) ® Xi (1)] 
from normalized orbitals 



(r) 



1 

L V 7T 

1 l~2 



L V 7T 



exp 



exp 



(x-X Q ) 2 + y 2 



{x + X ) 2 + y 



2 1 



L 2 



(Al) 

(A2) 
(A3) 



Here Xo-(l) and X<r(2) are spinors of the carriers 1 and 
2, and a =t,|, while r = (x,y). There are two vari- 
ational parameters: L and Xq. The former, width L, 
is introduced to lower the variational energy. The role 
of the latter, displacement Xq, is more important: the 
exchange interaction between holes and Mn spins takes 
place only for Xq ^ 0. This parameter introduces separa- 
tion of the holes with opposite spins, and is important for 
building the intuitive picture of effective interaction dis- 
cussed in the main text and below. However, the circular 
symmetry can be broken in any direction, for example, 
the x-axis can be the division line between Mn spins up 
and down, this would be reflected by the following re- 
placement 

(x±X ) 2 +y 2 ^x 2 + (y±Y ) 2 



in Eqs. (|A2IA3|) . 

The energy without the p - 



d exchange is 



-Enonm (L, Xq) 



2h 2 



*L 2 



-m*uj 2 (L 2 + 2X 2 ) (A4) 



Li 



X 



^Vtt Ry* ^ exp(-2^H I ( 2^ ) , 



X 



where L\ — v 2 Iq 



2h 



(A5) 



Here, m* is the effective mass, and Hluq is the energy 
quantum of the 2D harmonic oscillator. The second 
line of Eq. (|A4|) introduces the Coulomb repulsion of 
holes, integrated using 2D Fourier transforms [A[l[, Jo 
is the modified Bcsscl function of the 1st kind, Ry* = 



m*e 4 /[32 (ireh) 2 ] is the effective Rydberg, e is the elec- 
tron charge and e is the dielectric constant. We calculate 
the exchange energy using 



A / \ Ps (r)\d\ 



(A6) 



(Eq. (A7) in Ref. [42|). Here, A = x Mn N |/3| 5/2 is 
the exchange splitting, expressed through the fraction of 
cations replaced by Mn atoms, iMn, and the exchange 
constant N (3 < 0. The 2D spin density, p s , is expressed 
in h units. It is related to (sk) in the main text through 
Ps( r fc) /^z = ( s k) /3, where bar denotes spatial averaging 
over a cell Nk, and h z is the QD height. In the case of 
PS, the hole spin density is 



Pps 



[vl(r)-V d (r)] 12. 



(A7) 



One finds J ppsd 2 



owing to normalization of 



<-Pu,d (?*)• We obtain the exchange energy 



E ex (L,X ) =-A erf 



where erf is the error function. The results in Fig. (2a 
main text) are obtained by minimizing 



•Etot = E n 



{L,X ) + E ey , (L,X ) 



(A8) 



with respect to Xq and L. 

For iMn = 0, and without Colomb interaction, the 
wavefunction in Eq. (|A1[) with Xq = is the exact singlet 
groundstate. The wavefunction can be generalized to the 
case of elliptical-confinement (uj x < oj y ) by replacing L 
with the lengths L x>y = y/^lx^y, the latter are defined in 
the main text. The energy of the singlet is Huj x + huj y . 

The variational wavefunction for the triplet in a cir- 
cular QD is the same as in Eq. (4) in Ref. [AQ. The 
energy for the lowest triplet (T) in an elliptical QD 
with XMn = is obtained by placing one of the holes 
in the lowest excited orbital: p x . The energy is then 
Et = 2fvjj x + hoj y . Thus, the singlet-triplet gap decreases 
as huj x = hujo (1 + d) 2 with increasing d. The smaller 
this gap, the lower iM n value is sufficient for the triplet 
to become the groundstate (E ex for the triplet is always 
more negative than E cx for PS). This fact can be used 
to understand the dependence of the PS-T boundary on 
increasing d, Fig. (3 main text). 

Spin and charge densities for the triplet are propor- 
tional. Hence, the double-peaked spin density in Fig. (3b 
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main text) can be explained in the variational model as 
a direct consequence of the form of triplet wavefunction: 
One hole is in the s orbital, the other in the p x orbital. 



Al.l. Effective repulsion for Pseudosinglet 

It is instructive to fix L and consider the trends of 
particular contributions to i^tot [Eq. (|A8[) ] as a func- 
tion of Xq. Kinetic energy docs not depend on Xq, 
Eki n = 2h 2 / (m*L 2 ). Both Coulomb repulsion and Mn- 
cxchange terms decrease with increasing Xq , to zero and 
— Ao, respectively. The latter limiting value can be un- 
derstood as the energy of two independent magnetic po- 
larons with opposite spin projections (hole spin-up and 
Mn spins-down for x > 0, while hole spin down and Mn 
spins up for x < 0). This can be seen from pps (r) in 
Eq. (|A7j) . which has a minimum and a maximum sepa- 
rated, in this limit, by 2Xq, see also Fig. (3a main text). 
Confinement (potential) energy increases as m*cj 2 AQ, 
preventing the two holes from the complete separation, so 
that Xo < +oo, and l-E^x (L, X < oo)| < A . Thus, we 
can treat -dE cx (L, X ) /dX cx ±A exp [-2X1/ L 2 ] l L 
as an effective repulsive force between the two holes. This 
force enhances the Coulomb repulsion. 

The above reasoning remains qualitatively correct at 
non-zero temperatures, as long as they are are lower than 
the phase transition temperature. For T > 0, we calcu- 
late the variational exchange energy from a generalization 
of Eq. (|A6l) 



Ao / Ps (r) B s 



5 Psjr) 
2 h z k B T 



r. 



(A9) 



resulting from Eqs. (2, 6 main text). Eq. (|A6j) above 
is obtained in the limit of T — > 0, when the Brillouin 
function B s [|/3| 5p s (r) / (2h z k B T)} -> sign [p a (r)]. 

Before numerically exemplifying the above model of 
effective interaction, we notice that the numerical hole- 
spin density (s^) would not be a convenient choice for 
studying the effective repulsion in PS, as it vanishes in 
the small iMn limit. However, the spatial separation of 
the two holes can also be seen in the charge density (in 
units of |e|) 

g(r) = ip u (rf + tp d (r) 2 . 

The distance of any of the two symmetrical maxima from 
the QD center, r — 0, is given by 



-Rpeak = Xq tanh 



4i?pcakATo 



L 2 



in the variational framework. To satisfy this condition 
for i? P cak > 0, the condition X > L/2 must hold. We 
neglect Coulomb interaction for simplicity for the rest of 
this subsection. For iMn -> Owe find L — > L 1: [Eq. (|A5[) ]. 
and Xq — > 0, so there is a lower-limit XMm below which 
g has only one extremum; a maximum at r = 0. Above 




X Mn (%) 

Figure Al. (color online) Variational results for positions of 
maxima of charge density g for PS at T = 0. Solid green line: 
Rpcak, dotted blue line: Xo, dashed black line: L/2. We use 
the standard QD parameters defined in the main text, except 
for the effective Rydberg Ry* = 0. 



this XMn, there are two maxima at r = (±i? poa k, 0) and 
a minimum at r = 0. The results of these considerations 
are shown in Fig. jAQ] 



A1.2. Effective attraction for Triplet 

In the case of the triplets discussed in the main text, 
the hole spin density px has a constant sign, so that 
J p^d 2 r = ±1. Thus, at zero temperature, the maximum 
energy gain occurs for unidirectional alignment of Mn 
spins, and no changes of the wavefunction are needed to 
attain the maximum, cf. Eq. (|A6|) . This suggests that the 
T = limit is not helpful in the case of the triplets, since 
at any T > 0, the local exchange energy, \/3\ p s (r) /h z , is 
smaller than k B T for Mn spins sufficiently remote from 
the center, so that some localization (change of wave- 
function) is beneficial. As expected, the magnetization 
patterns for the triplets remain unidirectional at any tem- 
perature. 



A1.3. Effective interaction in the Exact 
Diagonalization method 

Figure A[2] shows the numerical results obtained by the 
exact diagonalization (EXD) method described in the 
main text. The smallest Mn-content used is iMn = 0.1%. 
For this value, the number of Mn spins in the disk-shape 
volume of radius L\ and height h z is N = 3. Hence, 
for the very low contents iMn < 0.1%, the continuum 
approximation would no longer be justifiable. Typical 
values of iMn in experiments arc higher than 0.1%. 
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X M n (%) 



Figure A2. (color online) Distance, 7? pe ak, from the center of 
a QD to the position of one of two symmetrical maxima in the 
charge density g, calculated in EXD. iMn > 0.1% as explained 
in the text here. Red: PS, black: triplet. Dashed: T = 0, 
solid T — 2 K. The oscillations at high iMn are artefacts due 
to numerical precision. We use the standard QD parameters 
defined in the main text, taking into account the Coulomb 
interaction. 

Same as in the variational approach, we find that in 
EXD the distance i? pC ak depends on iMn, justifying the 
notion of effective interaction. Note the opposite behav- 
ior of PS vs. triplet at T > 0, representing repulsion and 
attraction respectively, as discussed in the main text. 

Comparing the EXD results for PS with the variational 
ones (Fig. A[T]), we see that the off-center maxima start 
to appear at a much lower value of % n (0.6 vs. 1.6%). 
This is because the Coulomb interaction, present in EXD, 
favors separation of holes. 

Finally, we note that a different quantum approach to 
interaction of magnetic bipolarons was used in Ref. [AJU . 

A1.4. Core-halo pattern 

As explained in the main text, a different Mn-spin pat- 
tern forms for some Mn-doping profiles. This pattern can 
be studied variationally by replacing the tp u and (fd func- 
tions in Eq. (|Alj) by two Gaussians with Xq = 0, but with 
different widths: L u for spin up and for spin down, 
see Eq. (3) in Ref. [A[|[. 

A2. PHOTOLUMINESCENCE IN TYPE-II QDs 



holes are confined. Thus, the exchange interaction of Mn 
with electrons is much weaker than with holes. 




CB 



1 1 VB 

ZnSe (Zn,Mn)Te ZnSe 

Figure A3, (color online) Band alignment for a type-II band 
profile, found in, e.g., ZnSe/ZnTe QDs. CB and VB stand for 
conduction and valence bands respectively. 

Figure A[5] shows the transitions calculated in the 
manuscript. We assume that owing to fast energy re- 
laxation, luminescence from groundstatcs dominates over 
other initial states. Two types of transitions occur, in 
which the QD occupancies change as: 2 — >• 1 and 1—^0. 
These two kinds of lines are separated by an energy deter- 
mined mainly by hole repulsion (charging energy) ~ 35 
meV. In general, the final state of a 2 — > 1 is not the 
initial state of the subsequent 1 — > transition, since re- 
laxation of a Mn-spin pattern may occur between the two 
radiative events. The numerical calculations of the final 
states of the 2 — > 1 transitions are not selfconsistent, but 
rather carried out for Mn-magnetization rrik fixed in the 
pattern established when two holes are still present. 

As argued in the main text, interband transitions 
from p-like states of electrons may arise due to break- 
ing of the circular-symmetry by the Mn-pattern (the 
non-magnetic confinement potential remaining circularly 
symmetrical). The symmetry selection rules for the in- 
terband transitions are given by overlaps of electron and 
hole wavefunctions. Expanding the shifted Gaussian or- 
bitals, Eqs. (|A2IA3|) for small X /L, we obtain 

ip u (r) ~ s (r) + ^j-p x (r) , 

ip d (r) ~ s (r) - -j-p x (r) , 

where s = -\/2/7rcxp (— ? ,2 /L 2 ) /L, and p x — 
y/8/n' xcxp (-r 2 /L 2 ) / L 2 . Thus, for x M n > 0, the hole 
wavefunction contains p x orbitals, which have the same 
symmetry as electron p-likc orbitals. The hole wavefunc- 
tion for the core-halo pattern does not have the p orbitals 
in the equivalent expansion. 



In a (Zn,Mn)Te/ZnSe QD system, holes are confined 
in the QDs, while electrons are in the barrier (weakly 
bound to holes through Coulomb interaction) [A HQ. 
This is an example of the type-II band profile (Fig. A[3]), 
as opposed to the type-I profile, where both electrons and 
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B=0 B>0 




vac-*— t 

Figure A4. (color online) Jabloriski diagram [A0] of the tran- 
sitions discussed in the main text. Horizontal black lines: 
single and two-hole levels, as well as the zero-hole ground- 
state, "vacuum". Blue arrows: radiative transitions. Zig-zag 
arrows: non-radiative relaxation of Mn spins. Hollow arrows: 
spins of holes. An electron spin (not shown) must be an- 
tiparallel to hole spin to recombine [A @] . The energies of 
the radiative transitions in the left /right panel are shown in 
Figs. 4/5 of main text. Eps/Et is the pseudo-singlet/triplet 
energy. E s or E p are the energies of a single hole in an s-like 
or p-like state respectively. The number —1 is the angular 
momentum of the triplet and p-like states. 
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